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The many-parametric fit of the LEP2 data on e + e~ — ► e + e~, t + t~ processes is performed 

to estimate signals of the Abelian Z'-boson beyond the standard model. The model-independent 
relations between the Z' couplings to the standard model particles allow to describe the Z' effects 
in lepton processes by 4 independent parameters. No signal is found by the complete LEP2 data 
set, and the 1.3a signal is detected by the fit of the backward bins. The Z' couplings to the vector 
and axial-vector lepton currents are constrained. The comparisons with the one-parameter fits and 
with the LEP1 experiments are performed. 
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PACS numbers: 14.70.Pw, 13.66.-a 



I. INTRODUCTION 

The precision test of the standard model (SM) at LEP 
gave a possibility not only to determine all the parame- 
ters and particle masses at the level of radiative correc- 
tions but also afforded an opportunity for searching for 
signals of new heavy particles beyond the energy scale of 
it. On the base of the LEP2 experiments the low bounds 
on parameters of various models extending the SM have 
been estimated and the scale of new physics was obtained 
III 0] Although no new particles were discovered, a 
general believe is that the energy scale of new physics 
to be of order 1 TeV, that may serve as a guide for ex- 
periments at the LHC. In this situation, any information 
about new heavy particles obtained on the base of the 
present day data is desirable and important. 

A lot of extended models includes a so-called Z' gauge 
boson - a massive neutral vector particle associated with 
an extra U(l) subgroup of the underlying group. Search- 
ing for this particle in either model-dependent or model- 
independent approaches is widely discussed in the litera- 
ture Q. In the papers of the present authors 0, 0, 
the new method and the observables for the model- 
independent search for the Z'-boson were proposed and 
applied to analyze the LEP2 experiment data. In con- 
trast to other model- independent searches, our approach 
gives a possibility to pick out uniquely this virtual state 
and determine its characteristics. It is based on two 
ingredients: 1) The relations between the low-energy 
parameters are derived from the renormalization group 
(RG) equation for the fermion scattering amplitude and 
therefore called the RG relations. Due to them the num- 
ber of unknown Z' parameters to be measured consid- 
erably decreases. 2) The kinematic properties of the Z' 
signals in various scattering process exhibit themselves 
when these relations are accounted for. In these papers 
the one-parametric observables were introduced and the 
signals of the Z' have been determined at the la con- 
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fidence level (CL) in the e + e~ — > process, and 

at the 2(7 CL in the Bhabha process. The Z' mass was 
estimated to be 1—1.2 TeV that gives a good chance to 
discover the particle at the LHC. 

Recently the final data of the LEP collaborations DEL- 
PHI and OPAL 0,0 were published and new more pre- 
cise estimates could be obtained. In the present paper 
we update the results of the one-parameter fit and per- 
form the complete many-parametric fit of the LEP2 data 
to estimate a possible signal of the Z'-boson. Usually, 
in a many-parametric fit the uncertainty of the result 
incerases drastically because of extra parameters. On 
the contrary, in our approach due to the RG relations 
between the low-energy couplings there are only 2-3 in- 
dependent parameters for the LEP scattering processes. 
Therefore, we believe that an inevitable increase of confi- 
dence areas (CA) in the many-parametric space could be 
compensated by accounting for all the accessible experi- 
mental information. As it will be shown, the uncertainty 
of the many-parametric fit can be comparable with the 
uncertainty of the previous one-parametric fits in Refs. 
0, Q- In this approach the combined data fit for all 
lepton processes is also possible. 



II. THE ABELIAN Z' BOSON AT LOW 
ENERGIES 

Let us adduce a necessary information about the 
Abelian Z'-boson. This particle is predicted by a num- 
ber of grand unification models. Among them the Eq 
and 50(10) based models [9] (for instance, LR, x ~ ^ 
and so on) are often discussed in the literature. In all 
the models, the Abelian Z'-boson is described by a low- 
energy U(X) gauge subgroup originated in some symme- 
try breaking pattern. 

At low energies, the Z'-boson can manifest itself by 
means of the couplings to the SM fermions and scalars as 
a virtual intermediate state. Moreover, the Z-boson cou- 
plings are also modified due to a Z-Z' mixing. In prin- 
ciple, arbitrary effective Z' interactions to the SM fields 
could be considered at low energies. However, the cou- 
plings of non-renormalizable types have to be suppressed 
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by heavy mass scales because of decoupling. Therefore, 
significant signals beyond the SM can be inspired by the 
couplings of renormalizable types. Such couplings can be 
derived by adding the new U (l)-terms to the electroweak 
covariant derivatives D cw in the Lagrangian |4| 
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where 4> is the SM scalar doublet; Jl, /r are the SM 
left-handed fermion doublets and right-handed fermion 
singlets; Z M denotes the U(l) symmetry eigenstate; and 
2/0, yf L and j// B mean the unknown couplings characteriz- 
ing the model beyond the SM. Instead of the couplings to 
the left-handed and right-handed fermion states it is con- 
venient to introduce the couplings to the axial- vector and 
vector currents: a } = (y^ - y/J/2, v f = (y ft + 2//J/2. 

The spontaneous breaking of the electroweak symme- 
try leads to the Z-Z' mixing. In case of the Abelian 
Z'-boson, the Z-Z' mixing angle 9q is determined by the 
coupling j/0 as follows |5j 
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where 9\y is the SM Weinberg angle, and a cm is the elec- 
tromagnetic fine structure constant. Although the mix- 
ing angle is a small quantity of order m^, 2 , it contributes 
to the Z-boson exchange amplitude and cannot be ne- 
glected at the LEP energies. 

The Lagrangian leads to the following interactions 
between the fermions and the Z and Z' mass eigenstates: 
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are 



where / is an arbitrary SM fermion state; v 
the SM couplings of the Z-boson. 

In a particular model the couplings vt and af take 
some specific values. In case when the model is unknown, 
these parameters and the mixing angle remain potentially 
arbitrary numbers. However, this is not the case if one 
assumes that the underlying extended model is a renor- 
malizable one. As it was shown in Ref. some of 
them have to be correlated due to renormalizability. The 
corresponding relations are 



v f - a f = v f * - a f * 



(4) 



where /* is the SU(2) partner of a fermion /, and T^j 
is the third component of the fermion isospin. They are 
motivated by the renormalization group equations at the 



Z' decoupling energies and also connected with the U(l) 
gauge symmetry of the Lagrangian. These relations cover 
all the popular models of the Abelian Z' boson allowing 
the model-independent searches for this particle. 

The relations ||2J incorporate the most common fea- 
tures of the Abelian Z'-boson. As it is seen, the axial- 
vector coupling is universal for all the fermion flavors. 
So, in what follows we will use the shorthand notation 
a = a e = = a T . The axial- vector coupling determines 
also the coupling to the scalar doublet and, consequently, 
the mixing angle. As a result, the number of indepen- 
dent couplings is significantly reduced. Considering the 
lcptonic processes e + e~ — > £ + £~ (£ = e, /x, t), one has to 
keep 4 unknown couplings: a, v e , t> M , and v T . Moreover, 
the RG relations serve to uniquely specify a kinematic 
domain of deviations from the SM predictions due to the 
virtual Z' boson. Thereof a single definition of the Z' 
signal can be done. 

In our analysis, as the SM values of the cross-sections 
we use the quantities calculated by the LEP2 collabora- 
tions 

0HE1II3- They account for either the one-loop 
radiative corrections or initial and final state radiation 
effects (together with the event selection rules, which are 
specific for each experiment). The deviation from the 
SM is computed in the improved Born approximation. 
This accuracy is sufficient for our analysis, leading to the 
systematic error of fit results less than 5-10%. 

The differential cross-section of the process e + e~ — > 
£ + £~ deviates from the SM value by various quadratic 
combinations of couplings a, v e , v^, v T . For the Bhabha 
process it reads 



da 

dz 



da SM 
dz 



2 2 



where the factors are known functions of the center-of- 
mass energy and the cosine of the electron scattering an- 
gle z. The deviation of the cross-section of e + e~ — ► /t + /x - 
(t + t~) processes has the similar form 



da 

dz 



do 



SM 



dz 



(6) 



This is our definition of the Z' signal. 

Since the Z' couplings enter the cross-section together 
with the inverse Z' mass, it is convenient to introduce 
the dimensionless couplings 



mz 



-af, Vf 



Airmz 



■Vf, 



(7) 



which can be constrained by experiments. 

Note again, that the cross-sections in Eqs. ©-© ac- 
count for the relations Q through the functions fi(z), 
fz(z), fi(z), since the coupling y^ (the mixing angle 6>o) 
is substituted by the axial coupling constant a. Usually, 
this dependence on the scalar field coupling is neglected 
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at all, when a four-fermion effective Lagrangian is ap- 
plied to describe physics beyond the SM p|. However, in 
our case, when we are interested in searching for signals 
of the Z'-boson on the base of the effective low-energy 
Lagrangian (JIJ, these contributions to the cross-section 
are essential. 



III. MANY-PARAMETER FITS 

As the basic observable to fit the LEP2 experiment 
data on the Bhabha process we propose the differential 
cross-section 



da 



Bhabha 
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dz 



dz 
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where i runs over the bins at various center-of-mass en- 
ergies y/s. The final differential cross-sections measured 
by the ALEPH (130-183 GeVjnj), DELPHI (189-207 
GeV, H), L3 (183-189 GeV, [nf), and OPAL (130-207 
GeV, 0) collaborations are taken (299 bins). 

As the observables for e + e~ — » t + t~ processes, 

we consider the total cross-section and the forward- 
backward asymmetry 
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where i runs over 12 center-of-mass energies y/s from 130 
to 207 GeV. We consider the combined LEP2 data Q 
for these observables (24 data entries for each process). 
These data are more precise as corresponding differen- 
tial cross-sections. Our analysis is based on the fact that 
the kinematics of the s-channel processes is rather sim- 
ple and a differential cross-section is effectively the two- 
parametric function of the scattering angle. The total 
cross-section and the forward-backward asymmetry in- 
corporate a complete information about the kinematics 
of the process and therefore are an adequate alternative 
for the differential cross-sections. 

The data are analysed by means of the x 2 fit- Denoting 
the observables (JSJ) © by cr.;, one can construct the x 2 - 
function, 
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Sat 



(10) 



where <t cx and da are the experimental values and uncer- 
tainties of the observables, and a th are their theoretical 
expressions presented in Eqs. I©-©. The sum in Eq. 
(|10H refers to either the data on one specific process or 
the combined data on several processes. By minimizing 
the x 2_ function, the maximal-likelihood estimate for the 
Z' couplings can be derived. The x 2 -function is also used 
to plot the CA in the space of a, v e , v^, and v T . 

For all the considered processes, the theoretic predic- 
tions af 1 are linear combinations of products of two Z' 



couplings 



af 1 — CjjAj, 



(11) 
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where CV, are known numbers. In what follows we use the 



matrix notation a th — crj , <7 0X 
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Cij , A 



The uncertainties Sai can be substituted by a covariancc 
matrix D. The diagonal elements of D are experimental 
errors squared, Da — ((5of x ) 2 , whereas the non-diagonal 
elements are responsible for the possible correlations of 
observables. The % 2 -function can be rewritten as 



X 2 (A) 



,th\T 



D^(a c 
CA) T D-\a c 



a th ) 
-CA), 



(12) 



where upperscript T denotes the matrix transposition. 



The x -function has a minimum, xt 



at 
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corresponding to the maximum-likelihood values of Z' 
couplings. From Eqs. (|12fl . (|13|) we obtain 

X 2 (A)-xLn = (A-AfD-^A-A), 

D = (C^D^C)- 1 . (14) 

Usually, the experimental values cr ox are normal- 
distributed quantities with the mean values <r th and the 
covariance matrix D. The quantities A, being the su- 
perposition of er cx , also have the same distribution. It 
is easy to show that A has the mean values A and the 
covariance matrix D. 

The inverse matrix I) -1 is symmetric and can be di- 
agonalized. The number of non-zero eigenvalues is deter- 
mined by the rank (denoted M) of D~ l . The rank M 
equals to the number of linear-independent terms in the 
observables a th . So, the right-hand-side of Eq. (jTlfl is 
a quantity distributed as x 2 with M degrees of freedom 
(d.o.f.). Since this random value is independent of A, the 
CA in the parameter space (a, v e , w M , v T ) corresponding 
to the probability (3 can be defined as |l2| : 



X' 



^min 



+ XclA M )i 



(15) 



where Xcl is the (3- level of the x 2 -distribution with 

M d.o.f. 

In the Bhabha process, the Z' effects are determined by 
3 linear- independent contributions coming from a 2 , w 2 , 
and av e (M = 3). As for the e + e~ — > ,t + t~ pro- 

cesses, the observables depend on 4 linear-independent 
terms for each process: a 2 , v e Vfj_, v e a, av^ for e + e~ — > 
fi + /i~ ; and a 2 , v e v T , v e d, dv T for e + e~ — > t + t _ (Af = 4). 
Note that a number of terms in the observables for dif- 
ferent processes are the same. Therefore, the number 
of x 2 d.o.f. in the combined fits is less than the sum 
of d.o.f. for separate processes. Hence, the predictive 
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FIG. 1: The 95% CL areas in the (a, v e ) plane for the Bhabha, 
e + e~ — > )i + fi~ , and e + e~ — > t + t~ processes. 



power of the larger set of data is not drastically spoiled 
by increased number of d.o.f. In fact, combining the data 
of the Bhabha and e + e _ — > n + n~ ( T+T ~) processes to- 
gether we have to treat 5 linear-independent terms. The 
complete data set for all the lepton processes is ruled by 
7 d.o.f. As a consequence, the combination of the data 
for all the lepton processes is possible. 

The parametric space of couplings (a, v e , w M , v T ) is 
four-dimensional. However, for the Bhabha process it is 
reduced to the plane (a, v e ), and to the three-dimensional 
volumes (a, v e , v^), (a, v e , v T ) for the e + e~ — > /i + ^t~ and 
e + e _ — > t + t~ processes, correspondingly. The predic- 
tive power of data is distributed not uniformly over the 
parameters. The parameters a and v e are present in all 
the considered processes and appear to be significantly 
constrained. The couplings or v T enter when the pro- 
cesses e + e~ — > or e + e~ — > r + r _ are accounted 
for. So, in these processes, we also study the projection 
of the CA onto the plane (a, v e ). 

The beginning of the parametric space, a = v e = 0, 
corresponds to the absence of the Z' signal. This is the 
SM value of the observables. This point could occur in- 
side or outside of the CA at a fixed CL. When it lays out 
of the CA, this means the distinct signal of the Abelian 
Z' . Then the signal probability can be defined as the 
probability that the data agree with the Abelian Z' bo- 
son existence and exclude the SM value. This probability 
corresponds to the most stringent CL (the largest Xcl) a ^ 
which the point a — v e — is excluded. If the SM value 
is inside the CA, the Z 1 boson is indistinguishable from 
the SM. In this case, upper bounds on the Z 1 couplings 
can be determined. 

The 95% CL areas in the (a, v e ) plane for the separate 
processes are plotted in Fig. ^ As it is seen, the Bhabha 
process constrains both the axial-vector and vector cou- 
plings. As for the e + e~ — > and e + e~ — > r + r~ 
processes, the axial-vector coupling is significantly con- 
strained, only. The CAs include the SM point at the 
meaningful CLs, so the experiment could not pick out 
clearly the Abelian Z' signal from the SM. An important 
conclusion from these plots is that the experiment signifi- 
cantly constrains only the couplings entering sign-definite 
terms in the cross-sections. 
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FIG. 2: The projection of the 95% CL area onto the (a, v e ) 
plane for the combination of the Bhabha, e + e~ —> and 
e + e~ — » t + t~ processes. 

The combination of all the lepton processes is pre- 
sented in Fig. |21 There is no visible signal beyond the 
SM. The couplings to the vector and axial-vector elec- 
tron currents are constrained by the many-parameter fit 
as \v e \ < 0.012, \a\ < 0.018 at the 95% CL. If the charge 
corresponding to the Z' interactions is assumed to be of 
order of the electromagnetic one, then the Z' mass should 
be greater than 0.67 TeV. For the charge of order of the 
SM SU{2) L coupling constant m z > > 1.4 TeV. One can 
see that the constraint is not too severe to exclude the 
Z' searches at the LHC. 

Let us compare the obtained results with the one- 
parameter fits in Ref. 0. Fitting the current data with 
the one-parameter observable, we find the updated values 
of the Z' coupling to the electron vector current together 
with their la uncertainties: 



ALEPH 


vl = 


-0.11 ±6.53 x 10~ 4 


DELPHI 


vl = 


1.60 ± 1.46 x 10~ 4 


L3 


vl = 


5.42 ±3.72 x 10~ 4 


OPAL 


vl = 


2.42 ± 1.27 x 10~ 4 


Combined 


vl = 


2.24 ±0.92 x 10~ 4 . 



As it is seen, the most precise data of DELPHI and OPAL 
collaborations are resulted in the Abelian Z 1 signals at 
one and two standard deviation level, correspondingly. 
The combined value shows the 2a signal, which corre- 
sponds to 0.006 < \v e \ < 0.020. 

On the other hand, our many-parameter fit con- 
strains the Z' coupling to the electron vector current 
as \v e \ < 0.012 with no evident signal. Why does the 
one-parameter fit of the Bhabha process show the 2a CL 
signal whereas there is no signal in the two-parameter 
one? Our one-parameter observable accounts mainly 
for the backward bins. This is in accordance with the 
kinematic features of the process. Therefore, the dif- 
ference of the results can be inspired by the data sets 
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FIG. 3: The 68% CL area in the (a, v e ) plane from the back- 
ward bins of the Bhabha process in the LEP2 experiments 
(the shaded area). The hatched area is the 68% CL area from 
the LEP 1 data on the Bhabha process. 

used. To check this, we perform the many-parameter 
fit with the 113 backward bins (z < 0), only. The \ 2 
minimum, Xmin = 103.0, is found in the non-zero point 
\a\ = 0.00017, v e = 0.015. This value of the Z' cou- 
pling v e is in an excellent agreement with the mean value 
obtained in the one-parameter fit. The 68% CA in the 
(a, v e ) plane is plotted in Fig. There is a visible signal 
of the Abelian Z' boson. The zero point a — v e — (the 
absence of the Z' boson) corresponds to x 2 = 107.9. It is 
covered by the CA with 1.33<r CL. Thus, the backward 
bins show the 1.33cr signal of the Abelian Z' boson in 
the many-parameter fit. So, the many-parameter fit is 
less precise than the analysis of the one-parameter ob- 
servables. 

At LEP1 experiments E3 the Z-boson coupling con- 
stants to the vector and axial- vector lepton currents (<?y, 
<7a) were precisely measured. The Bhabha process shows 
the ler deviation from the SM values for Higgs boson 
masses m H > 114 GeV (see Fig. 7.3 of Ref. jig). This 
deviation could be considered as the effect of the Z-Z' 
mixing. It is interesting to estimate the bounds on the 
Z' couplings following from these experiments. 

Due to the RG relations, the Z-Z' mixing angle is com- 
pletely determined by the axial- vector coupling a. So, the 
deviations of gv , 9a from their SM values are governed 
by the couplings a and v e , 

9v - g v M = -49.06aw e , g A - gf 1 = 49.06a 2 . (16) 

Let us assume that a total deviation of theory from ex- 
periments follows due to the Z-Z' mixing. This gives an 
upper bound on the Z' couplings. In this way one can 
estimate is the Z' boson excluded by the experiments or 
not. 

The la CL area for the Bhabha process from Ref. E3l 
is converted to the (a,v e ) plane in Fig. The SM val- 
ues of the couplings correspond to the top quark mass 
m t — 178 GeV and the Higgs scalar mass tuh = 114 
GeV. As it is seen, the LEP1 data on the Bhabha process 
is compatible with the Abelian Z' existence at the lcr CL. 



The axial-vector coupling is constrained as \a\ < 0.005. 
This bound corresponds to a 2 < 2.5 x 10~ 5 , which agrees 
with our previous one-parameter fits of the LEP2 data for 
e+e- -> fi+fi- , t+t- processes @ (a 2 = 1.3±3.89 x 10~ 5 
at 68% CL). On the other hand, the vector coupling con- 
stant v e is practically unconstrained by the LEP1 exper- 
iments. 



IV. DISCUSSION 

LEP collaborations have reported about a good agree- 
ment between the experimental data and the predictions 
of the SM 0, H EH 113 0ur analysis of the leptonic 
processes based on the same data set and the same SM 
values of cross-sections showed that the existence of Z' 
boson with the mass of order 1-1.2 TeV is not excluded 
at the l-2cr CL. We observed this in one-parameter fits 
0, Q and in the many-parameter fits in the present pa- 
per. The estimated Z' parameters derived in different 
methods are in good agreement with each other. So, we 
are faced with a necessity to find a possible explanations 
of this discrepancy. We believe that the reason is in the 
RG relations, which play a crucial role in treating exper- 
imental data. As we showed, the RG relations served to 
reduce the number of unknown parameters and extract a 
maximal information about the signal of the particle of 
interest from the experimental data set. If the RG rela- 
tions are not taken into account, no signals of Z' will be 
found. 

LEP collaborations performed also model-dependent 
fits concerning popular Z' models. Since these models 
are included in the ones suiting the RG relations (0J, it 
is interesting to compare their analysis with our results. 
In experiments Q, H [I EH E3 the low bound on the Z' 
mass was obtained. It varies from 400 to 800 GeV at the 
95% CL dependently on the Z' model. These bounds 
allow the Z' boson with the mass of order 1 TeV, being 
compatible with our results. A possibility to select Z' 
signals in specific scattering processes was not discussed 
in these papers. This is the reason why no signals were 
observed. 

In our analysis we treat the data for leptonic processes 
only. LEP2 collaborations measured also the total cross- 
sections of the electron-positron annihilation into quark- 
antiquark pairs. The Abelian Z' signal in e + e~ — > qq 
process is characterized by 5 independent parameters (for 
example, two Z' couplings to electron, a and v e , and 
three Z' couplings to d, s and b quarks). There are 8 
linear independent terms in the cross-section. As one can 
check, the experimental statistics of 12 cross-sections for 
different center-of-mass energies is completely insufficient 
to constrain significantly the parameters of the Z' boson. 

As we have shown in Ref. pj, there is the 2a signal of 
the Abelian Z' boson in the one-parameter fit of LEP2 
data for the Bhabha process. This result is reproduced 
in the present paper by fitting the updated experimental 
data. In the present analysis we applied many-parameter 
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fits for leptonic processes for different sets of bins in- 
cluded. In particular, for the backward bins (responsible 
for the signal due to the kinematics of the process) the 
1.3cr signal of the particle is found. The fit of the com- 
plete set of bins constrains the Z' couplings to vector and 
axial- vector electron currents allowing the Z 1 boson with 
the mass of order 1 TeV. Thus, we have to conclude that 
the LEP2 data allow the existence of quite light Z'-boson, 
which has a chance to be discovered in the nearest future, 
at the LHC. We believe that the RG relations used in the 
present analysis will be also important in searches for the 
Z' boson at the LHC. 
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